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Abstract —A major performance degrading factor in free space 
optical communication (FSO) systems is atmospheric turbulence. 
Spatial diversity techniques provide a promising approach to 
mitigate turbulence-induced fading. In this paper, we study 
the error rate performance of FSO links with spatial diversity 
over atmospheric turbulence channels described by the Double 
Generalized Gamma distribution which is a new generic sta¬ 
tistical model covering all turbulence conditions. We assume 
intensity modulation/direct detection with on-off keying and 
present the BER performance of single-input multiple-output 
(SIMO), multiple-input single-output (MISO) and multiple-input 
multiple-output (MIMO) FSO systems over this new channel 
model. 

Index Terms —Atmospheric turbulence. Double GG distribu¬ 
tion, bit error rate, error rate performance analysis, Free-space 
optical systems, spatial diversity, MIMO. 

I. Introduction 

F REE-SPACE optical (ESO) communication enables wire¬ 
less connectivity through atmosphere using laser trans¬ 
mitters at infrared bands. These systems provide high data 
rates comparable to fiber optics while they offer much more 
flexibility in installation and deployment. Since they operate 
in unregulated spectrum, no licensing fee is required leading 
to a cost-effective solution Q-El. 

A major performance limiting factor in FSO systems is 
atmospheric turbulence-induced fading (also called as scintil¬ 
lation) 141. Inhomogenities in the temperature and the pressure 
of the atmosphere result in variations of the refractive index 
and cause atmospheric turbulence. This manifests itself as 
random fluctuations in the received signal. In the literature, 
several statistical models have been proposed to model this 
random phenomenon. Historically, log-normal distribution has 
been the most widely used statistical model for the random 
irradiance experienced over atmospheric channels iH-iiQi. 
This model is however restricted to weak turbulence conditions 
and has large deviations from the experimental data when the 


strength of turbulence increases. 

In an effort to come up with a more general model to cover 
a wide range of turbulence conditions, other statistical models 
have been proposed in the literature which include the K 
ifTH . I-K ifT^ . log-normal Rician ifTSll . Gamma-Gamma lfT4ll . 
M II15II and Double Weibull m distributions. In our recent 
work, we proposed the so-called Double Generalized Gamma 
(Double GG) as a unifying distribution for the irradiance 
fluctuations im. This model is valid under all range of 
turbulence conditions and contains most of the existing models 
in the literature as special cases. 

In our previous work ini, as an initial performance study, 
we derived the BER performance of single-input single¬ 
output (SISO) FSO link over Double GG channels. In this 
work, we extend our performance analysis to multiple-input 
multiple-output (MIMO) FSO systems. MIMO FSO systems 
are known to mitigate turbulence-induced fading and sig¬ 
nificantly improve the performance. Some earlier results on 
MIMO FSO systems over log-normal, K, negative exponential 
and Gamma-Gamma channels can be found in ISl-Ga. In 
this paper, we study the error rate performance of single¬ 
input multiple-output (SIMO), multiple-input single-output 
(MISO) and MIMO FSO systems employing intensity mod¬ 
ulation/direct detection (IM/DD) with on-off keying (OOK) 
over independent and not necessarily identically distributed 
(i.n.i.d.) Double GG turbulence channels. 

The rest of the paper is organized as follows: In Section II, 
we introduce the MIMO FSO system model. In Section III, 
we provide the BER expressions for SIMO, MISO and MIMO 
FSO links. In Section IV, we present numerical results to con¬ 
firm the accuracy of the derived expressions and demonstrate 
the advantages of employing spatial diversity over SISO links. 
Finally, Section V concludes the paper. 
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II. System Model 


derivation) 


We consider an FSO system employing IM/DD with OOK 
where the information signal is transmitted via M apertures 
and received by N apertures over the Double GG channel. 
The received signal at the n* receive aperture is then given 
by 
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III. BER Performance 


The optimum decision metric for OOK is given by 
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Therefore, the average error rate can be expressed as 
/ 




MIMO 


where x represents the information bits and can be either 0 or 
1, is the Additive White Gaussian noise (AWGN) term with 
zero mean and variance = Nq/2 , and 77 is the optical-to- 
electrical conversion coefficient. Here, Imn is the normalized 
irradiance from the m* transmitter to the n* receiver whose 
pdf follows ins Eq. (1)] 
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is the Meijers G-function defined in 
Eq.(9.301)], p and q are positive integer numbers that satisfy 
p/q = 71/72 and A{j;x) A x/j ,..., {x P j - l)/j , and 
Pi > 0.5 is a shaping parameter modeling the severity of 
fading. The distribution parameters 7 ^ and VLi, i = 1, 2 , of 
the Double GG model can be identied using the following 
equations 


where fi (I) is the joint pdf of vector I = (/n, / 12 , •. •, Imn)- 
The factor M in ® ensures that the total transmitted powers of 
diversity system and SISO link are equal for a fair comparison. 
On the other hand, the factor N is used to ensure that sum 
of the N receive aperture areas is the same as the area of 
the receive aperture of the SISO link. The integral expressed 
in ([H does not yield a closed-form solution even for simpler 
turbulence distributions; however, it can be calculated through 
numerical multi-dimensional integration. Similarly, we can 
use multidimensional Gaussian quadrature rule (GQR) 
techniques to calculate the BER for MISO case, i.e.. 
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which does not yield a closed form expression either. 

In the following, we focus on the SIMO case and investi¬ 
gate the BER performance under the assumption of optimal 
combining (OC) with perfect CSI where the variance of the 
noise in each receiver is given by cr^ = No/2N. Therefore, 
replacing M = 1 in ® we obtain 
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where and ay are respectively normalized variances of 
small and large scale irradiance uctuations. 

The cumulative distribution function (cdf) of Double GG 
distribution can be derived from (|2|i as ifTTll 


Eq. ([Tot does not yield a closed-form solution and requires 
N-dimensional integration. Nevertheless, the Q-function can 
be well-approximated as Q{x) ~ e~A ~/12 -|- /4 

and thus the average BER can be obtained as 
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where r = (ri, r 2 ,..., rjv) is the received signal vector. 
Eollowing the same approach as 113, Eol, the conditional 
bit error probabilities are given by (see 03 for details of 


The above integral can be evaluated by first expressing the ex¬ 
ponential function in terms of the Meijer G-function presented 
in lIZTl eq. (11)] as 

exp(-x) = Go;° [t lo ] ( 12 ) 

Then, a closed-form expression for (fTTT l is obtained using lIZTl 
Eq. (21)] as 
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where A (n, u) is defined in (fT4l i at the top of the next page. 
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The derived BER expression in (fTsT i for SIMO FSO system 
with OC can be seen as a generalization of BER results over 
other atmospheric turbulence models. Eor example, if we insert 
7 i = 1 and 17^ = 1 in (fO]) . we obtain the BER expression over 
Gamma-Gamma channel. Setting /3i = 1 in (fOl l. we obtain 
the BER for Double Weibull channel. On the other hand, for 
7 i = 1, Oi = 1 and (32 = 1, (fOl l reduces to (12) of lfT9l 
reported for the K-channel. Appendix provides the details on 
these. 

As an alternative to OC, we also consider equal gain com¬ 
bining (EGC) where the receiver adds the receiver branches. 
In this case, the average BER can be expressed as 


PSIM0,ECG = J fl (I) Q £ /n j dl (18) 

It should be noted that (fTsT l is equivalent to ® obtained for the 
MISO ESO links. Another method is selection combining (SC) 
which is the least complicated of the combining schemes since 
it only processes one of the diversity apertures. Specically, the 
SC chooses the aperture with the maximum received irradiance 
(or electrical SNR). Therefore, the pdf of the output of SC 
receiver can be obtained as 
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The average BER can be then expressed as 
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which can be efficiently calculated through numerical means. 



Fig. 1. Comparison of the average BER bet\veen SISO and different diversity 
techniques for plane wave assuming i.i.d. turbulent channel defined as channel 
b. 



Fig. 2. Comparison of the average BER between SISO and different diversity 
techniques for plane wave assuming i.i.d. turbulent channel defined as channel 
c. 


IV. Numerical Results 

In this section, we present BER performance results of 
SIMO, MISO and MIMO ESO systems over Double GG 
channels and quantify the performance improvements over 
SISO systems. We consider the following four scenarios of 
atmospheric turbulence conditions reported in ini. 





























• Channel a: Plane wave and moderate irradiance fluc¬ 
tuations with 7 i = 2.1690, 72 = 0.8530, mi = 0.55, 
m 2 = 2.35, 111 = 1.5793, ll 2 = 0.9671, p = 28 and 
q=ll 

• Channel h: Plane wave and strong irradiance fluctu¬ 
ations with 7 i = 1.8621, 72 = 0.7638, mi = 0.5, 
m 2 = 1.8, Hi = 1.5074, H 2 = 0.9280, p = 17 and 
9 = 7. 

• Channel c; Spherical wave and moderate irradiance 
fluctuations with 71 = 0.9135, 72 = 1.4385, mi = 2.65, 
m 2 = 0.85, Hi = 0.9836 and H 2 = 1.1745, p = 7 and 
q = 11. 

• Channel d: Spherical wave and strong irradiance fluc¬ 
tuations with 7 i = 0.4205, 72 = 0.6643, mi = 3.2, 
m 2 = 2.8, Hi = 0.8336 and H 2 = 0.9224, p = 7 and 
( 7 = 11 . 

Figs. 1-2 present the average BER over i.i.d. channels 
defined as channel b and channel c respectively. For SIMO 
FSO links employing OC receivers, we present approximate 
analytical results which have been obtained through (fOl l along 
with the Monte-Carlo simulation of (fTOl) . As clearly seen from 
Figs. 1-2, our approximate expressions provide an excellent 
match to simulation results. As a benchmark, the average 
BER of SISO FSO link is also included in these figures. It is 
observed that multiple receive aperture deployment signicantly 
improves the performance. For instance, at a target bit error 
rate of 10“®, we observe performance improvements of 47.2 
dB and 67.2 dB for SIMO FSO links with N = 2 and 3 receive 
apertures employing OC with respect to the SISO transmission 
over channel b. Similarly, for channel c, at a BER of 10“®, 
impressive performance improvements of 51.5 dB and 64.3 
dB are achieved for SIMO links with N = 2 and 3 employing 
OC compared to the SISO deployment. It is also illustrated 
that EGC receivers yield nearly the same performance as OC 
receivers. For example, in Fig 2, for N = 2 the performance 
difference between OC and EGC receivers is merely 0.4 dB 
at a BER of 10“^. Also as expected, EGC and OC receivers 
outperform SC counterpart. 

Figs. 3-4 demonstrate the BER performance of SIMO 
FSO links employing OC, EGC and SC receivers over non- 
identically distributed (i.n.i.d.) Double GG channels. Similar 
to i.i.d. results, our approximate closed-form expressions again 
yield nearly identical match to simulation results. We further 
compare the performance of i.n.i.d. case with respect to i.i.d. 
case presented in Figs 1-2. For example, to achieve a BER of 
10“® in SIMO links with N = 2 over i.n.i.d. channels a and b, 
we need 8.2 dB, 8.5 dB and 4.9 dB less in comparison to i.i.d. 
case respectively for OC, EGC and SC receivers. Note that in 
Fig. 1, we assume that both of the two channels between the 
transmitter and receivers are described by channel b. Thus, 
since the channel a is less severe than the channel b, we need 
less SNR in comparison to i.i.d case to obtain the same BER. 
On the other hand, to achieve a BER of 10“® for SIMO links 
with N = 2 over i.n.i.d. channels c and d, we need 6.8 dB 
more for OC and EGC receivers and 11.6 dB more for SC 



Fig. 3. Comparison of the OC, EGC and SC receivers for SIMO FSO links 
over two i.n.i.d. atmospheric turbulence channels defined as channel a and 
channel b. 



Fig. 4. Comparison of the OC, EGC and SC receivers for SIMO FSO links 
over two i.n.i.d. atmospheric turbulence channels defined as channel c and 
channel d. 

receiver in comparison to i.i.d channels. Note that in Fig. 2, 
both of the two channels between the transmitter and receivers 
are described by channel c. Therefore, as the channel d is more 
severe than the channel c, we need more SNR in comparison 
to i.i.d case to achieve the same performance. 

V. Conclusions 

In this paper, we have investigated the BER performance of 
FSO links with spatial diversity over atmospheric turbulence 
channels described by the Double GG distribution. We have 
obtained an efficient and unified closed-form expression for 
the BER of SIMO FSO systems with OC receiver which 

















generalizes existing results as special cases. For MISO and 
MIMO systems, we have presented BER performance based 
on numerical calculations of the integral expressions. Our 
numerical results have demonstrated that spatial diversity 
schemes can significantly improve the system performance 
and bring impressive performance gains over SISO systems. 
Our comparisons among SIMO FSO links employing OC, 
EGC and SC receivers have further demonstrated that EGC 
scheme presents a favorable trade-off between complexity and 
performance. 


Appendix 

Proposition 1: Inserting 7 ^ = 1, = 1 and (3i = 1 in ( 1731 ). 

the BER of SIMO FSO links using optimal combining over 
K-channel is obtained. 

Proof: First by replacing 7 ^ = = p = <7 = /3i = 1 in (IThl l 
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which coincides with (21) of ifT^ . 

Proposition 2: Inserting 7 ^ = 1, = 1 in ( ITil ). the BER 

of SIMO FSO links using optimal combining over Gamma- 
Gamma channel is obtained. 

Proof: First by replacing 7 ^ = = p = q = 1 in (fTSI l and 

dnll, we obtain = l/T (/3i,„) F (/32,„) and = I3f)^l3f^^. 
Then by inserting all the values as well as / = 1 and k = 2 
in (HI), we obtain 
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Proposition 3: Inserting [3i = 1 in f ITil ), the BER of SIMO FSO 
links using optimal combining over Double-Weibull channel is 
obtained. 

Proof: Inserting /3i = lin (fThl l and (fTTl i, we obtain 
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